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Abstract 
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interest, such as vector nonlinear Schrodinger-type equations and three reso- 
nant wave equations are considered. 
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1 Introduction 



The theory of solitons originated long time ago with the discovery of the integrabil- 
ity of the Korteweg-de Vries (KdV) equation and with the method of the spectral 
transform to solve it pQ. Many other soliton equations, with various degrees of 
novel mathematical features and applicative interests (see, for instance, [2j and jHj), 
were then found by extending and generalizing the method of the spectral transform. 
Among these integrable equations, the most notable one is the nonlinear Schrodinger 
(NLS) equation (subscripted variables denote partial differentiation and x and t are 
the independent variables): 



where 7 is the real dispersion parameter. The great importance of this equation 
is due to its integrability and universality [3] which are related consequences of 
a perturbative multiscale analysis of (a large class of) dispersive nonlinear wave 
equations (see also 0). 

The key property of integrable evolution Partial Differential Equations (PDEs) is 
that they express the condition that two linear differential equations (Lax pair), 
both for the same unknown function and whose coefficients depend also on a com- 
plex (spectral) parameter, are compatible with each other. The compatibility of 
these two equations is a local condition and it provides local properties of the asso- 
ciated integrable PDE, in particular local conservation laws (continuity equations) 
and Hamiltonian structures. The construction of solutions of integrable nonlinear 
evolution PDEs is, however, a different matter. In fact, in addition to the initial 
value, one has to specify the domain of the space variable x together with the val- 
ues which the solution (and/or its x-derivatives) takes on the domain boundary. 
Once the (appropriate) boundary values are fixed, then one should make use of the 
Lax pair to solve the Cauchy initial value problem or, more modestly, to construct 
special solutions. In both cases, various methods of solution have been devised de- 
pending on the boundary conditions imposed on the solution, and several problems, 
for certain initial-boundary values, remain open. 

In this paper we give the explicit construction of special solutions based on Darboux- 
Dressing Transformations by means of algebraic and local operations only, with no 
need to referring to specific boundary values. The construction of soliton solutions, 
with their appropriate boundary conditions, will be considered in a subsequent paper 
for the case in which the space domain is the entire real axis, —00 < x < +00. 
Historically, the initial-value problem for the KdV equation in this space domain 
with vanishing values at the boundary, x = ±00, was the first to be solved. Then, 
by the same spectral technique but applied to a different Lax pair, Zakharov and 
Shabat were able to solve the initial-value problem for the NLS equation (^) on 
the entire real x-axis and for both vanishing = 6] and nonvanishing [7] values at the 
boundary x = ±00. The possibility to solve these two distinctive cases turned out 
to be quite relevant in nonlinear optics where the soliton behaviour in both cases 





s = 



±1, 



(1) 
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has been observed. If the boundary values are vanishing, soliton solutions of (JIJ 
exist for s = —1, while if the boundary values are nonvanishing (and appropriately 
prescribed) stable soliton solutions of exist for s = +1, while for s = —1 solitons 
exist but are unstable against small perturbations Solitons are usually referred 
to as "bright" solitons in the first case (vanishing boundary values) as they are 
light-pulses in a dark background, and as "dark" solitons in the second case (non 
vanishing boundary values) as they are pulses of darkness in a light background. In 
the latter case a further distinction is made between "black" and "grey" solitons. 
For a guide to the vast literature on solitons see, for instance, 0. 

The propagation of pulses in nonlinear media is more generally modelled by multi- 
component fields, and this requires to generalize the NLS equation (JTJ) to a system 
of coupled nonlinear Schrodinger equations. The coupling between the field compo- 
nents, which may represent different polarization amplitudes or fields with different 
(resonating) frequencies, as dictated by physical contexts, may or may not lead to 
integrable nonlinear equations. Though physical models may have terms which make 
the propagation equation nonintegrable, still the investigation of sufficiently close 
integrable equations, because of our good analytical control of them, is certainly 
worth and valuable also in an applicative context. Integrable systems of coupled 
NLS equations, termed vector nonlinear Schrodinger (VNLS) equations as the de- 
pendent variable u(x, t) is a D-dimensional vector, u = {u^ l \ ■ ■ ■ , w^), are known 
with various coupling terms. For instance, for D = 2 the system 

uP =t 1 [u^-2s(\u^\ 2 + \u^\ 2 )u^], s = ±l, U 

was first introduced by Manakov [TOj- Similarly to the scalar NLS equation (£[]), 
also this system possesses bright-solitons (i.e. with vanishing boundary values) for 
s — — 1 and (stable) dark-solitons (i.e. with nonovanishing boundary values) for 
s = +1. This system (j2J) remains integrable even if generalized to 

«« = i 7 [ w « - 2(sM 1) \ 2 + s 2 \u^\ 2 )u«] , , , 

i4 2) =i 7 [«S- 2(s 1 \u^\ 2 + s 2 \u^\ 2 )u^}, Sj = ±l, J = l,2; 

in fact, only three systems have been found integrable for two NLS equations coupled 
by square modulo terms of this form, namely equation (jSJ) with s\ = s 2 = ±1 and 

Sl = s 2 = i nn. 

Different ways of coupling two or more Schrodinger equations in a nonlinear way 
which preserves integrability have been discovered. Some of them are included in the 
class of wave equations considered here and are displayed in section 2. Additional 
systems of wave equations which do not posses the Schrodinger-type second order 
dispersion term since they are only first order in the ^-derivatives, and therefore 
they may or may not model dispersive waves, are also treated here. The most 
notable example of such systems is the one describing the coupling of three waves 
at resonance (see section^, a model which applies to several physical contexts such 
as optics, fluid dynamics and plasma physics. 
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All systems of coupled Partial Differential Equations (PDEs) discussed in this paper 
are particular reductions of a general matrix evolution equation, namely the matrix 
partial differential equation (j3J) . The integrability of this "mother" equation, which is 
guaranteed by the associated Lax pair, entails the possibility to solve the initial value 
problem by the spectral transform (or inverse scattering) method. However, the task 
of setting up the formalism and finding out the relevant integral equation has not 
been carried out completely. Indeed, this task is certainly easier if one assumes that 
the solution vanishes at the boundary, while in the case of nonvanishing boundary 
values it is rather involved because of the many branch points which occur in the 
complex plane of the spectral variable. Even in the case of the .D-dimensional 
VNLS (j21j) . which is a simple reduction of the system (@J), the spectral method 
in the case of nonvanishing boundary values, well understood for the scalar NLS 
equation Q [T2|, has not been taken to a state which is amenable to easy access in 
applications (the case D = 2 has been recently considered in and partial results 
have been obtained in [14 J. On the other hand, in the VNLS case, the dark vector 
soliton phenomenology has actracted considerable attention among nonlinear optics 
scholars, and explicit soliton expressions have been derived by direct methods (for 
instance Hirota method) JE1> rather than by using the Lax pair. 

The aim of this paper is to provide a general method of explicit construction of 
soliton solutions of the mother equation (|1J), and therefore of all the various wave 
propagation equations which obtain from it by reductions. This method makes use 
of the Lax pair and goes via the Darboux-Dressing Transformations (DDT) |19j-[2*2^ 
in a standard way; see also . However, the nonvanishing of the boundary values 
at infinity in the variable x introduces novel features with respect to the more usual 
context of vanishing boundary values. In fact, the novelty here regards more the 
case of dark solitons, rather than that of bright solitons. 

In the next section |2] we introduce the general formalism and display some of the 
reduced equations which we deem representative of potential applicative models. 
Some of these wave equations follow from a generalization which has been recently 
introduced in [21] and with the purpose of modelling solitons which behave as 
boomerons (or trappons) (see also [2H] and [27]). Section El is devoted to the DDT 
technique in its general setting. Section 13. 1[ deals with DDTs characterized by a 
complex (non real) pole, while section I3~2l deals with DDTs with a real pole, these 
two cases being appropriate to cope with different boundary values. In the last 
section 0] we add comments and remarks. 
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2 General formalism and multicomponent wave 
equations 



All systems of coupled PDEs considered in this paper are special (reduced) cases of 
the following matrix PDE: 



3^ 



Qt = [C(°\Q}+a[C^,Q x ]-a{Q,W}-t 1 a(Q xx -2Q ,. 
W X = {CM,Q*}, W 

where the dependent variables Q = Q(x,t) and W = W(x,t) are (N^ + ' + N^) x 
_|_ n(-)} block matrices of the form 

o=f 0j %r„ Q<+> )• w=( n wM °«->V (5) 



In self-evident notation, the diagonal entries W /< - + - ) , are square matrices of 

dimension x and N^~' x while the two off-diagonal rectangular 

blocks Q( + \ respectively Q^~\ are x N^~\ respectively x N^ + \ matrices, 
and being arbitrary, positive integers. Here, and in the following, the 
constant diagonal matrix a in (j3J) is, in self-evident notation, 

a _ | ljV(+)xAr(+) OjV(+)xAT(-) I ^ 
\Qn(-)xN(+) -IjVC-)xAT(-) / 

while C^ -*, C^ 1 -* are arbitrary constant block-diagonal matrices, 

00) = V« r «'-'), ,- = 0,1, (7) 

and the block-diagonal matrix W = W(x,t) is an auxiliary dependent variable. 
As usual [A, B] and {A, B} are the commutator AB — BA and, respectively, the 
anticommutator AB + BA. The constant coefficient 7 is the real dispersion pa- 
rameter, 7 = 7*. The systems of coupled NLS equations which are obtained 
by reduction of equation (J3J) have been considered in the special case with 
7 = having been extensively reported in [23] and j2H|, where it is shown that 
their single soliton solutions feature boomeronic and trapponic behaviours. It is 
worth noticing that setting = ca, and therefore W = 0, the linear terms 
[C(°\Q] and a[C^\ Q x ] can be both transformed away by the obvious transforma- 
tion Q(x,t) -> Q(x,t) = exp(-C (0) t)Q(x - 2c£, t) exp(C (0) t) and one obtains the 
matrix evolution equation 

Qt = -ijv(Q xx - 2Q 3 ) , (8) 

that is, the standard matrix version of the NLS type equation. However, whenever 
is a generic full block-diagonal matrix (see (|7ji). equation (j3J) is the most general 
second order differential equation which genuinely generalizes the standard equation 
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(JBJ). Indeed, it allows to treat the somewhat simpler class of evolution equations 
corresponding to the choice 7 = 0, the canonical form of which only features first 
derivatives with respect to the space variable x (rather than second derivatives as is 
the case of NLS type equations). The simplest system of equations of this type can 
in fact be reformulated so as to coincide with the standard equations describing the 
resonant interaction of 3 waves (see below) with the remarkable consequence that 
its solutions may feature the boomeronic or trapponic behavior, as recently reported 
in and £?]. 



2.1 Lax Pair and reductions 



The matrix equation (J1J) is the compatibility condition for the Lax pair 

ip x = Xiff, iff t = Tip, (9) 
where ip, X and T are (iV (+) + iV (_) ) x (JV« + JVH) 

square matrices, ip = ip(x, t, k) 
being a common solution of the two linear ordinary differential matrix equations © 
while X = X(x, t, k) and T = T(x, t, k) depend on the coordinate x , the time t and 
the complex spectral parameter k according to the definitions 

X(x,t,k) = -ikJ + Q(x,t), (10a) 



T(x, t, k) = -2jk [-ikJ + Q(x, t)} - 2ikC^ + i^a [Q 2 (x, t) - Q x (x, £)] + 

(10b) 

-aW + a[C^\Q] +C<°>. 
In these expressions the diagonal traceless constant matrix J is 

iV (+) - 2 /JVH \ 

N (+) + N n' 1 N (+) + N n'yo _iv(+) J ■ ( U} 

In order to simplify the notation, from now on (and as already done in the previous 
formula) we do not specify the dimension of the matrices and 1, and we may even 
omit to write the matrix 1 altogether, as we trust the reader will not be confused 
by this omission. 

Let us now consider the condition 

Q\x,t) = SQ{x,t)S (12) 

on the solution Q(x, t) of the matrix evolution equation Q where the dagger stands 
for hermitian conjugation. The constant matrix 5* is block diagonal, 

s =(tV->); (13) 
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its off-diagonal blocks are vanishing rectangular matrices while its diagonal blocks 
S« and are, respectively, x N^ + ' and ' x ' diagonal matrices 

whose diagonal elements Sn, with no loss of generality, are signs, namely 



5 (±) = diag (s\ 



(±) 



,(±) 
'jv(±) 



1 . 



(14) 



This of course implies the relations S 2 = 1, S^ 2 = 1, ^ 



The reduction equation (fT2^) is well motivated by the fact that it captures several 
interesting models of dispersive propagation of multicomponent waves in weakly 
nonlinear media. In order to support this claim, we display some of such model 
equations in the next subsection. To this aim it is convenient to rewrite the matrix 
PDEs (gj) in terms of the blocks Q ( ~\ W i+) and W { ~\ see ©. These read 



Q(±) =C r(0)(±)Q(±)_Q(±) C -(0)(= F ) ± 



"Cf(i)(±) gW-g^cWW 

2g(±)gwg(±) 



{±) 

^xx 



(15a) 



W K±) = [cW(±), Q(±)gW]. (15b) 

where respectively are the JVW x N^ + \ respectively x TV^, 

constant square matrix blocks of see (|7|). The reduction condition (|12J1 is 

accounted for by introducing the dependent variable U (x, t) through the definitions 

Ql-\x,t) = U(x,t), Q i+) {x,t) = S i+) U\x,t)S { - ) . (16) 

Since the constant 7 is real, these expressions of and in terms of the single 
variable U(x,t) are compatible with the equations (fT5Jl which then reduce to 



u t = c(°x-) c/ -ucw - [cwh U x - U x C««l + 

+ [W<-> U + U W< + >] + i7 [C/ ra -2U W U] , 



(17a) 



(17b) 



fH=[C«H,[7S(+)[7t5H]. (17c) 

Here C/ is an iV*") x iV( + ) rectangular matrix, (see ()16j) and (J3J)), while the auxiliary 
variables t) are square matrices, respectively is an iV( + ) x matrix 

and is an iV' ' x ]\f' ' matrix, and it is easily seen that they satisfy the 

"hermitianity" conditions 

W {+) = _ S (+) w (+)\ S {+) ; W {-) = _ S (-) w {-)\ $(-) ; (18) 

that is 

W\x,t) = -SW(x,t)S. (19) 
Similarly, the constant matrices C^)^ satisfy the following conditions: 

CW(+) = _(-y5W^(+)t £(+), j = 0,1, (20a) 
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C^'>H = - S H C {j) ^ , j = 0, 1. (20b) 

As for the "sign" matrices and S^~\ see ([14)1 . we note that one could set, for 
instance, = 1 with no loss of generality, but we prefer to keep the symmetrical, 
though redundant, notation (JHj). 

2.2 Multicomponent wave equations: cases of interest 

We observe that the general matrix nonlinear evolution equation (j!7|) may well 
specialize itself to quite a large family of coupled wave equations by playing with 
various choices of the integers and iV* - ** , and of the constant matrix coefficients 
C^'^ + ' and C^~\ This exercise is certainly worth doing since it happens that 
among this family of wave equations there are some which look interesting models 
in different applicative contexts. With this purpose in mind, and before proceeding 
further on our general setting, we briefly list here a few examples of model PDEs 
which we deem of interest in applications (for further details see also and |25j). 

Let us begin with choosing N^ +S) = 1 and N^' = D, C^'^' = and let us set S^ + ' = 
1 and S(~> = S =diag(si, • • • , sp) with = 1. In this case the dependent variable 
U(x,t) is a (column) D-dimensional vector u(x,t) = (u^'(x,t), • • • ,u^ D \x,t)) and 
the resulting VNLS equation 

= ij(vi xx - 2 < u, Su > u) , (21) 

is a simple generalization of the D = 2 case ©• Here and hereafter we adopt the 
bracket notation to indicate the (nonsymmetrical) scalar product of two vectors, 
namely 

D 

<u,v>=^n ( " ) V" ) . (22) 

n=l 

This system of NLS equations can be further generalized by adding coupling terms 
which originate from non vanishing coefficients C^'^K A simple instance of such 
generalization obtains for = 1 and = 2, together with = 1 and 

=dia K ( S\, S2), where s\ 2 — 1- Here the choice of the matrix coefficients is 

C (0)( + ) = C (i)( + ) = Qj C (0)(-) = ^ ^ s 1 a ^ ^ = ^ -6 ^ ^ (23) 

while the dependent variables are the two components (u^(x,t) , u^(x,t)) of the 
2-dimensional (column) vector U and the function w(x,t) is defined through the 
off-diagonal matrix 



u, 



D 



it 



(n)|2 



U 
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while = 0. With these specifications the VNLS equations read 

u t (1) = Sl auP> + bu x l) - si W 2 ) + 27 \u ( Ji - 2 U \u^\ 2 + s 2 \u^\ 2 ^) uP> 
u {2) = -s 2 a*uW - bu x 2) + s 2 W 1 ) +i 7 [ui 2 2 - 2 (si |u«| 2 + s 2 \u^u^ 

w x = 2 b Si s 2 m (1) w (2) * , 

(25) 

where a is an arbitrary complex coefficient and b is an arbitrary real coefficient. 
We note that this system results from a mixing of Schrodinger-type dispersion and 
quadratic nonlinearity as it occurs in the 3-wave resonant interaction. This exten- 
sion is simple and it might look trivial; however, because of the non commutative 
character of the matrix coefficients C^-°^~' and , which entails the introduc- 

tion of the new auxiliary dependent variable w(x,t), it introduces the interesting 
phenomenology of boomerons, i.e. solitons which have different asymptotic veloci- 
ties at t = +oo and t — — oo. These special soliton solutions were first introduced 
long time ago in [2E| (see also [2H]) in connection with the hierarchy of the ma- 
trix Korteweg-de Vries equation, but only recently boomerons appeared again in 
geometry [3U| and optics [27], [21]. 

Let us consider next = 2, = D, C^W = 0, =diag(si,s 2 ) and 

SH = S = diagf si, • • • , sp). In this case the matrix U has two columns, namely 
two D-dimensional vectors, u and v. They satisfy the following coupled VNLS 
equations 

u t = i"/(u xx - 2s[ +) < u, S^u > u - 2s 2 +) < v, S^u > v) , , , 

v t = ij{v xx - 2s[ +) < u, S^v > u - 2s { 2 +) < v, SH V > v) . ' " ' 

This system may of course be rewritten as one equation for a 2D-dimensional vector. 
However, this vector equation would be different from the equation which obtains 
in the case = 1, = 2D. Of course, similarly to the previous case, also this 
system can be generalized so as to feature boomeronic-type effects. 

Consider now the subcase D = 2, namely = N^~' = 2. The matrix U is 

2x2 square and therefore it can be given a different representation by using Pauli 
matrices: 

U = m (4) 1 + u {l) a 1 + u {2) a 2 + u {3) a 3 ; (27) 

with this parametrization of U the four dependent scalar variables are the com- 
ponents u^\j = 1, • • • ,4, of the 4-dimensional vector u = (u^ 1 ', u^ 2 \ u^'). 
However, the resulting system of evolution equations for these fields takes a vec- 
tor covariant form by specializing the "sign" matrices to = 1,S^ = si with 
s 2 = 1. In this case, the evolution equation for the 4-dimensional complex vector u 
has been first introduced in EU, and reads, in self-evident notation, 



Ut = ij{u xx — 2s [2 (u* ■ u)u — (u ■ u)u*]} . (28) 
Let us note that, via the transformation u^' — > vV' exp(i0j) with (pj four arbitrary 
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real constants, this system (}2e|) of four coupled NLS equations takes the form 

(29) 



4 fe) = %1 \ uf x - 2s 



vj=i / \i=i 

where the real constants = 5j k * = —Skj, with 5j k = 2(<f)j — 4>k), are three arbitrary 
phases. We also note that further reductions of this system are obtained by merely 
letting one, two or three components of the 4-vector u vanish. In fact, by setting 
for instance = 0, the system (|2U|) becomes a 3-dimensional VNLS equation 
with two arbitrary phases, say 5 2 i and 5 32 . Setting also = leads then to the 
2-dimensional case, which is 

= i 7 { u « - 2s[(\uW\ 2 + 2\u^\ 2 )u^ + exp(iS)u®W>]} , 
u f = %1 { u ^ - 2s[{2\u^\ 2 + \u^\ 2 )u^ + exp(-iS)uM 2 uW*]} , 

with only one arbitrary constant phase, 5 2 i — ^- This system has been first intro- 
duced in [32] with 5 = 0, while the two special cases with 5 = and 5 = n had 
been identified as the only integrable ones in the class of equations having the form 
written above, but with the two factors exp(i5) and exp(— iS) replaced by two apriori 
arbitrary real constants (US]- The case in which three components of u are vanish- 
ing yields the scalar NLS equation Larger systems of coupled nonlinear NLS 
equations are similarly obtained by starting with a bigger square matrix U(x,t). 
For instance 9-dimensional VNLS equations result from expanding the 3x3 matrix 
U in the basis obtained by adding the unit matrix to the SU(3) group generators. 
Other ways of coupling NLS equations may be deviced by different parametrizations 
of the rectangular matrix U. 

Finally, let us note that the well known 3-wave resonant interaction equation (in 
both its "explosive" and "non-explosive" versions) is a special case of our family of 
reduced equations as it merely coincides with the a = 7 = case of the system (|23|) 

(1) l. (1) (21 

u\ — bux = —Si w u [ ' , 

u f + bu x 2) = s 2 w*uW, (31) 

w x = 2bs 1 s 2 u^u {2 > . 
Moreover, having its physical applications in mind, one should also redefine for this 
equation the two independent variables x and t rather as "time" and, respectively, 
"space", with the implication that the three characteristic group velocities are 1/6, 
— 1/6 and 0. Indeed, this observation has played a basic role in the discovery of 
soliton solutions of the 3WRI equation such as boomerons, trappons, simultons [2Zj 
and pair creation. 



3 The Darboux-Dressing Transformation 

Let us now turn our attention to the method of construction of special solutions 
of the general system (JT7J). We first note that the reduction conditions (JTSJ), (fT9"j) . 
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(PD|). together with the expressions (jlOaj) and ()10b|) of the matrices X(x,t, k) and 
T(x, t, k) in the two linear equations (Lax pair) Q entail the following relations: 

X*{k*)Z + ZX{k) = 0, T f (F)£ + ET(fc) = , (32) 

with 

S = <^=(T-°-)) ' J? = 1 ' (33) 

As for the notation used in ()32|). we have omitted writing the dependence on the 
variables x and t and we maintain this omission in the following whenever it will 
cause no confusion. The property (}3*2*|) allows one to express the reduction condition 
induced by (fT2|) . (f"TT?|). (|2T?|) on the solution ip{k) of the two linear equations Q; this 
condition is given by the following equation: 

^(x, t, k*)Ei/;(x, t, k) = A(k } k*) , (34) 

where the matrix A(k, k*) is constant, namely x- and t-independent. Therefore, it 
is plain that the value of A(k,k*) depends only on the arbitrary value tp(xo,to, k) 
that the solution ip takes at a given point (xo,to) of the (x,t) plane. 

Consider now a second pair of matrices Q^(x,t) and W^(x,t), and assume that 
they have the same block structure of Q and W, see ©, and satisfy the same 
reduction conditions (fT2J) . (fK?J) . Let ip^(x,t, k) be a corresponding nonsingular (i.e. 
with nonvanishing determinant) matrix solution of 

V>£ 0) = X<V 0) , 4 0) = T(°V (0) (35) 

with X^ ' (x, t, k) and (x, t, fc) having the expressions (jlUa|) and (|10b|) with Q and 
replaced by Q and ll / (°' ) . Assume also that the initial condition i/j^°'(xq, t , k) 
is so chosen that the constant matrix A^(k, k*), where of course (see (jHH)) 

A^\k,k*) = ^\x,t,k*)^ \x,t,k), (36) 

coincides with A(k, k*), i.e. A^ \k, k*) = A(k,k*). Since both compatibility con- 
ditions, ip^) = ipte and ip xt = ijj tx , are satisfied, Q^(x,t), W^(x,t) and Q(x,t), 
W(x, t) are two different solutions of the same matrix evolution equations (jlj), and 
therefore it follows that the matrix 

D(x,t,k) = ifj(x,t } k)(ilj {0 \x,t,k))- 1 (37) 

satisfies the differential equations 

D x =XD- DX {Q) , D t = TD- DT {0) , (38) 

together with the algebraic (reduction) equation 

D\k*)ZD(k) = Z. (39) 
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The proof of these propositions is straight. 

The definition (|37|) can be viewed as a transformation of ip^ into ip 

ip(x, t, k) = D(x, t, k)^ (x, t, k) , (40) 

which consequently yields a transformation of and into, respectively, Q and 
W. Therefore, the dressing approach requires in the first place an explicit knowledge 
of Q(°\x,t), W^°'(x,t) and ip^°'(x,t,k). The next step is the construction of the 
transformation matrix D(x,t, k) via the integration of the ODEs (|38|). This task is 
however not straight since the coefficients X and T of these differential equations 
depend on the unknown matrices Q and W (see (llUaj) and (jlObjl ). The way of 
solving this problem goes through the a priori assignment of the dependence of the 
transformation matrix D(x,t, k) on the spectral variable k. 

In the following we will investigate the set of fc-dependent matrices D(k) which i) 
have a rational dependence on the complex variable k and ii) have nonvanishing 
k — > oo limit. Moreover, if we consider a rational dependence on k which can be 
factorized as product of simple-pole terms, we need to deal only with matrices D{k) 
which take the following one-pole expression 

D(x,t,k) = l + ^±, (41) 
k — a 

where the matrix R(x, t) is the residue at the pole k = a and the value of D(k) in the 
k —>■ oo limit is taken to be the identity for the sake of simplicity. The transformation 
(J4*U|) characterized by the matrix (j4Tj) has received considerable attention in the lit- 
erature JH] [22] ( see also (23] )• We refer to it as Darboux-Dressing Transformation 
(DDT), and its existence in our setting is proved below by construction. 

The way to obtain an explicit expression of the residue matrix R(x, t) depends on 
whether the pole a is off the real axis, a ^ a*, or on the real axis, a = a*. Therefore 
we treat these two cases separately. 



3.1 The Darboux-Dressing Transformation: complex pole 

Let us begin with considering the case in which a is not real, a ^ a*, while the 
other case, a real, is discussed in the next subsection 13.21 The starting point is 
the requirement that the matrix D(x,t, k), see (141 J) , satisfies the algebraic condition 
(I39j) and the differential equations ()38jl . The algebraic condition entails the two 
(equivalent) equations 

W+ ^ = , fl f E _^ = , (42) 
a — a* a — a* 

whose solution is 

R(x,t) = (a - a*)P{x,t) , (43) 
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where the matrix P(x, t) is a projector with the "hermit ianity" condition 

P 2 = p , pt = SPS _ (44) 

As for the differential equations replacing D(x,t,k) with its expression (|4*T|) 
(and yields the algebraic relations 

Q = gt°)+i(a-a*)[<7,P], (45a) 



W = W {0) -i(a-a*) , {a, P}} , (45b) 

which give the "dressed" matrices Q and W in terms of the "bare" matrices Q^°\ 
and the projector P, together with the two differential equations 

P x = X(a)P - PX {0) (a) , P t = T(a)P - PT (0) (a) , (46) 

whose integration goes as follows. Consider first the differential equation ()46)1 with 
respect to the variable x and replace Q with its expression ([45 aj) . The resulting 
equation is then the following nonlinear equation 

P x = X {0) (a)P - PX {0) (a) + i(a - a*)(JP - PJP) . (47) 

Let z be an eigenvector of P, and differentiate with respect to x the eigenvalue 
equation 

Pz = z . (48) 
By replacing then P x with the right-hand side of (|47j) . one arrives to the equation 

(1-P)[z x -X®(a*)z} = 0, (49) 

which implies that the vector z x — X^°'(a*)z is in the subspace on which P projects. 
At this point we may well assume that this subspace is one-dimensional. Indeed, it 
is easy to prove that, if P projects on a subspace of higher dimension n > 1, then 
the corresponding matrix D is a product of as many matrices of the form (j41j) 
with (gSD, 

D^,t,k) = l + {a ' a ' )P " >{X - t] , i=l,...,», (50) 

k — a 

as the dimension n of this subspace, all of course with the same pole in a, and 
all with Pw) projecting on a one-dimensional subspace. Therefore, with no loss of 
generality, we let P in the DDT matrix 

D(x,t,k) = l + ia -f )P ^ t \ (51) 

k — a 

project on the one-dimensional subspace of the vector z, with the implication, see 
()49p . that the vector z x — X^°'(a*)z is proportional to z. On the other hand, since 
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the vector z is defined here modulo a scalar factor function, we may choose this 
factor in such a way that z satisfies the differential equation 

z x = XM(a*)z. (52) 

The differential equation ()46|) with respect to the variable t can be treated in a 
similar way. The substitution of Q and W with their expressions ()45|) yields the 
nonlinear equation 

p t = T(°)(a)P - PT(°\a) + 2 7 (a - a*)(Q<®P - PQ®P)+ , , 

- 2i 1 (a 2 -a* 2 )(JP-PJP) + 2i(a-a*)\C ( < l \P}P. 1 j 

By differentiating now the eigenvalue equation (f4~S|) with respect to t and using both 
the equations (}4*5j) and (p)3j) . one ends up with the equation 

(1-P)[z t -T®(a*)z} = 0, (54) 

which, by the same arguments as above, implies that the vector z(x, t) satisfies the 
differential equation 

z t = T^(a*)z. (55) 

Once the two equations ()52|) and ()55|) are solved, the DDT transformation matrix 
D(x, t, k) is finally given by ([5Tf with 



this expression being implied by the algebraic conditions ()44|). At this point we 
conclude that the method of construction of a novel solution Q(x,t), W(x, t) of 
the evolution equations (0J), starting from the knowledge of given {seed) solution 
Q(°\x,t), W(°\x,t) is explicitly given by (j43|) with (|5^|) where the vector z(x, t) is 

z (x,t) = ?p (0) (x,t,a*)z . (57) 

Here ip^(x, t, a*) is the solution ip^°>(x,t,k) of the differential equations (jHB]) (Lax 
pair corresponding to Q^°\x, t), W^°'(x, £)), for k = a*, and it is assumed to be 
known, while z is an arbitrary constant (N^ + ^-dimensional vector. 

The construction of the novel solution Q(x,t), W(x, t) provides, via the formulas 
(jHJ) and the reduction (JTHJ), the construction of the novel solution U(x,t), W^(x, t) 
of the matrix evolution equation (|17|) which is our main concern here. The relevant 
expressions are obtained by first writing the block structure of the (iV( + ) + N^) x 
(]V r (+) _|_ j\K™)) projector matrix P(x,t) 

which therefore entails, in self-evident notation and by using (f43j) . the relations 

U = U i0) - 2i(a - a*)B , (59a) 
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W (+) = w (o)(+) -2i(a- a*)[CW + \ , (59b) 



W {-) = w (o)(-) + 2i( a - a*)[C (1)H ,5 H ] • (59c) 

In view of its use in future computations, we give to this formula a more explicit 
expression by using the form ("5~oj) of the projector P. To this aim it is convenient 
to split the vector z in two block column vectors, namely 



(60) 



where the vectors and z^ have dimension and, respectively, N^~> . Then, 
by inserting this block form of z in the diadic expression of P, (|56|) (and recalling 
(J!"""")), we arrive to the final relations 

U = um - 2i(a - a') < , w - _ < ^ m - , (61a) 



, w , ccww z(+)zMt£(+) - z (+) z (+)t 1 s , (+) r7( 1 )(+) , i 

w(+) = w(0)(+) _ 2if a-a*)^- - - - - - - - - (61b) 



3.2 The Darboux-Dressing Transformation: real pole 

Let us now investigate the alternative case in which the pole of the Darboux-Dressing 
matrix D(x, t, k), see (JHZJ), is real, a = a*. The way to treat this case is the same as 
that one we followed in the previous case, but the resulting equations to be solved are 
indeed different. Thus we first ask that D(x, t, k), as given by the general expression 
dlH) , satisfies both the algebraic condition ([59*]) and the differential equations (jUSJ). 
The algebraic condition implies two equations for the residue matrix R: 

ER + R^ = 0, R^R = 0. (62) 

These entail the following form of R 

R(x,t) = ip(x,t) P(x,t)E , (63) 

together with the conditions that the scalar function p(x, t) is real, the projector 
matrix P is hermitian, 

p = p * , p 2 = p , p = pt ; ( 6 4) 
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and it satisfies the equation 

PEP = . (65) 
Therefore, in the present case, the transformation matrix reads 

D(x, t) = l + ip(x, t) P fr ^ E . (66) 

k — a 

By the same arguments we have used in the previous case, one can show tha P may 
be assumed, with no loss of generality, to project on a one-dimensional subspace, 
namely 

zz^ , . 

P = — — , 67 

< z,z> 

where the vector z(x, t), because of the equation (ffiKJ) . is constrained by the orthog- 
onality condition 

<z,i:z>=0. (68) 

Consider now the differential equations ()38|) and insert in them the expression (}6T)j) . 
Since k is of course an arbitrary complex variable, we obtain the relations 

Q = Q {0) +p[P,<r]E, (69a) 



W = W {0) + p[C (1) ,{a,P}£], (69b) 

which give the novel solution Q(x,t), W(x,t) of (jH) in terms of the supposedly 
known solution Q^(x,t), W^°'(x, t), the function p(x,£) and the projector P(x,t). 
There also follows the relation 

T(a)-T (0) (a) = 4jap[a 7 P]E + 2z 7 p[Q (0) , PE] + 2i 7 p 2 yPE -2pa[C^\ aPS] , (70) 

which is not an independent relation but it is instrumental in deriving the formula 
(|72j) displayed below. In addition, one obtains two following differential equations, 
one with respect to x: 

( P P) X = P (x®(a)P + PX<M(a) + p 9 - f P] , (71) 



and one with respect to t: 



(pP)t = P (rW(a)P + PT(°)t(a) - 4a 7 p^P - 2i 1P - f P + 2py/) . (72) 

In these last three equations we have conveniently introduced the functions f(x,t), 
g(x,t), h(x,t) and m(x,t) according to the definitions 

f=<z,z>, g=<z,Sz>, h=< 5,EQ (0) z > , m=<l,SC (1) 2> . (73) 
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We have obtained the right-hand side of the differential equation (J72j) by using the 
relation (|7Uj). At this point we differentiate the eigenvalue equation (see (jHZj)) 

Pz = z. (74) 

When this is done with respect to x, one obtains the equation 

P[z x + X^\a)H] = z x - X^\a)z + - p^j z , (75) 

which implies that its right-hand side is proportional to the vector z, 

z x -X^{a)z+{^-p 9 -^z = pz. (76) 

On the other hand, the vector z, which has been introduced through the diadic 
expression (|67|). is defined only modulo a factor scalar function, and therefore, by 
taking advantage of this freedom, one can ask that 

P = P-t, 77 

P f 

with the implication that the vector z satisfies the differential equation 

z x = X®(a)z. (78) 
We also observe that the function p, via its own definition 

P[z x + X < - \ a yz}=pz , (79) 
takes the following expression 

11 f 

p=j<z, z x + X^( a y Z >=j(<z,z x > + <z x ,z>) = J -j , (80) 

where we have taken into account the equation (f7S|) and the definition of / (|75j). It 
is now readily found that combining this equation with the relation ()77|) entails that 
the function p(x, t) satisfies the following differential equation 



-g 



Let us now differentiate the eigenvalue equation (|71j) with respect to t. Following 
the same strategy as before, one concludes that the vector z satisfies the differential 
equation 

z t = T^(a)z, (82) 
while the following relations hold true 

P[z t + T^{a)H] = uz , (83) 
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v = EL + 2- (2a^g + i<yh - m) . (84) 
P f 

Again, the equation (JHHj) implies the expression 

v = j < z , z t + TW(a)+i >=j(<z,z t > + <z t ,z>) = ^, (85) 

which, together with ()84j) . yields the differential equation 

7 , 

' 1 2(2a-fg + i-fh-m) (86) 

for the function p(x, t) with respect to the variable t. Since the vector z(x, t) satisfies 
the two (compatible) differential equations (f7H|) and (|82|). its general expression is 

£(x,t) = ^ (0) (x,t,a) 5 , (87) 

where £o is an arbitrary constant (N^ + iV^^-dimensional vector. The two 
differential equations (JHTj) and (|86|h which are also compatible with each other (the 
proof is straightforward and it is not reported here), can be easily integrated since 
the functions f,g,h and m are known (see (|75jl ). The expression of their general 
solution then reads 

' SSSj + 2 /to dt ' [2«7^(^o, f ) + i-yh(x , t>) - m(x ,t')\ - /* dx'^rr', t) 
where Xo , to and p(xo,to) are arbitrary real constants. 

We conclude that, if Q^°'(x,t), W^°\x,t) and ip^°\x,t,k) are known, the explicit 
expressions of p(x,t) and z(x,t) given above yield, via and (JB7J) . the DDT 
matrix and therefore the new solution Q(x,t), W(x,t) through (|69|). As in the 
previous case, the corresponding formulas which give the expression of the solutions 
U(x, t), W^'(x, t) of the equation (fT7|) follow from the block structure of Q and W, 
of E, of the projector P 

B<+) fit 



and of the vector 5 



I ■ (90) 
The relevant relations then read 

U = U {0] + 2pBS i+) , (91a) 

W {+) = W (0)(+) + 2 p[C'( 1 )W, flW^W] , (91b) 



^(-) = W (0)(-) + 2 p[C7( 1 )(-), BH 5 H] ; 
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(91c) 



or, equivalently (see(|S7)l. (J89"j) and (JSDJ)) and more explicitly, 





(92a) 



< > + < z(-),|H > ' 




(92b) 




(92c) 



where the expression of p = p(x,t) is given by (|%5j) with (see (J7HJ)) 



f(x, t) =< |(+), |C+) > + < zH, > , 

=< fW, £(+)*(+) > + < tfHzH > , 
h(x,t) =< zM,U<MSl-W-) > - < f(-),5H[/(°)#) >, 
m(x,t) =< z(+),5(+)C , ( 1 )Wz(+) > - < zH^HcWH^-) > 



(93) 



We end this section noticing that the explicit formulas derived here and in the 
previous subsection 13.11 are meant to serve as the main tools to construct soliton- 
and, by repeated application of DDTs, multisoliton-solutions of the matrix equations 
(JT7J). However, these formulas have been obtained by algebra and local integration 
of differential equations. Therefore, the two important properties of the solutions 
they yield, namely their boundary values and their boundedness, are left out of our 
present discussion. These two issues, which are of basic relevance in applications, 
will be taken up in a following paper [33] where explicit expressions of solutions will 
be displayed and investigated. 



4 Conclusions and remarks 

In this paper, we have set up the formalism of the Darboux-Dressing Transform 
(DDT) to construct a new solution of the coupled matrix nonlinear integrable evo- 
lution PDEs 0, 



from a given known solution of the same equations. The choice of this system of 
equations, which has been first introduced in [21], is mainly motivated by the fact 
that it captures several interesting models of dispersive propagation of multicompo- 
nent waves in weakly nonlinear media. In fact, it should be pointed out that setting 
= ca implies that the linear terms [C^°',Q] and cr[C^ 1 ' ) , Q x ] can be both trans- 
formed away with the consequence that W = thereby obtaining the well known 
matrix NLS type equation 



Qt = [tf<°>, Q] + <t[CM,Q x ] - a{Q, W} - i 1( j{Q xx - 2Q 3 ) , 
W x = [CW,Q 2 ], 



Qt = -i^Qxx - 2Q 3 ) , 
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which originates well known VNLS equations. In contrast, whenever is a generic 
matrix, it is precisely the non commutative character of the matrix coefficients 
and which introduces the interesting phenomenology of boomerons, i.e. solitons 
which have different asymptotic velocities at t = +00 and t = — 00. These special 
soliton solutions were first introduced long time ago in connection with the hierarchy 
of the matrix Korteweg-de Vries equation, but recently boomerons appeared again 
in geometry and optics and are attracting a growing interest in applications. More- 
over and interestingly enough, the matrix evolution equation given above, which is 
the most general second order differential equation generated by the Lax pair Q, 
includes the simpler class of evolution equations corresponding to the choice 7 = 
which turns out to be nonlinear and non trivial. This class features first derivatives 
with respect to the space variable x (rather than second derivatives as is the case 
of NLS type equations), and the simplest system of equations of this type coincides 
with the standard equations describing the resonant interaction of 3 waves - yet their 
solutions feature, somewhat unexpectedly, the boomeronic or trapponic behavior. 
This is one instance of wave equations which are both in our class and of well known 
applicative relevance. We trust the interested reader will be attracted by looking for 
other wave equations via reduction of our "mother" system (@J) either by exploring 
the examples of such equations reported in or by engaging himself /her self in 
deriving new equations. 

The class of evolution equations we have chosen contains in fact, as reduced cases, 
systems of wave equations which, because of their dispersion (mainly of Schrodinger- 
type) and/or because of their nonlinearities, are, or promise to be, interesting models 
of applicative importance and it is with applications in mind that the DDT technique 
has been investigated and presented here almost as an algorithm to find explicit 
solutions. In its simplest version, this technique goes via a Darboux transformation 
of the matrix solution of the Lax equations, which introduces just one pole in the 
complex plane of the spectral variable. Though this method is known, we show that 
new algorithmic features arise if the pole is on the real axis. This case, on the other 
hand, is the one which matters in the construction of solutions with non vanishing 
boundary values (as for dark solitons). The relation between the reality of the pole 
and the boundary values, together with explicit solutions, will be discussed in a 
separate paper. 
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